OTIC  FILE  COPY  AD- A 146  175 


AD 


CRDC-CR*84087 


FLUID  MOTION  IN  A  ROTATING  AND 
NUTATING  CONTAINER 


by  THORWALD  HERBERT 

Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 

Blacksburg.  Virginia  24061 


JULY  1984 

US  Army  Armament,  Munitions  &  Chemical  Com  v<nd 
Aberdeen  Proving  Ground,  Maryland  21010 


»  0  9  28  048 


Disclaimer 


The  findings  in  this  report  are  not  to  be  construed  as  an 
official  Department  of  the  Army  position,  unless  so  designated  by 
other  authorized  documents.  , 


Disposition 


For  classified  documents,  follow  the  procedures  in  dod 
5200. 1-R,  Chapter  IX  or  or  DoD  5220. 22-M,  "Industrial  Security 
Manual,"  paragraph  19.  For  unclassified  documents,  destroy  by  any 
method  which  precludes  reconstruction  of  the  document. 


Distribution  Statement 

Approved  for  public  release;  distribution  unlimited. 


UNCLASSIFIED 

gC0WTrgaTOATi6N 


1*  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

2*.  SECURITY  CLASSIFICATION  AUTHORITY 

2b.  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 

4.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

CRDC-CR-84087 


REPORT  DOCUMENTATION  PAGE 

i  lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION  /AVAILABILITY  OF  REPORT 

Approved  for  public  release, 
distribution  unlimited. 

5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL 

Dept,  of  Engineering  Sci  &  Mech-  (if  applicable ) 

Institute1011^3  Pol^teclin:'-c 

6c  ADDRESS  (City.  Staff,  and  ZIP  Code) 

Blacksburg ,  VA  .24061 


7a.  NAME  OF  MONITORING  ORGANIZATION 


7b.  ADDRESS  (C/ty,  Staff,  and  ZIP  Code) 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Chemical  R&D  Center 


8b.  OFFICE  SYMBOL 
DRSMC-CLJ- 
IR  (A) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

DAAG29-81-D-0100  j •'.  /  ,  1 


8c  ADDRESS  (Gty,  State,  and  ZIP  Code) 

Aberdeen  Proving  Ground,  MD  21010-5423 


I  10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  PROJECT  H 

ELEMENT  NO.  NO.  I 


PROJECT 

TASK 

NO. 

NO. 

1L161102 

A71A 

WORK  UNIT 
ACCESSION  NO. 


It.  TITLE  (Include  Security  Classification) 

Fluid  Motion  in  a  Rotating  and  Nutating  Container 


12.  PERSONAL. AUTHOR(S) 

Herbert,  Thorwald 

13a.  TYPE  OF  REPORT 

Contractor _ 

16.  SUPPLEMENTARY  NOTATION 
COTR:  Miles  Millc 


13b.  TIME  COVERED  14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

I  frqmJu!  82  to  Jan  83|  198  4  July  52 


COTR:  Miles  Miller,  DRSMC-CLB-PA  (A)  (301)  671-2186 

17. _ _ COSATI  CODES _  18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

field _ group  sub-group  Projectile  stability 

15 _ 02 _ Liquid  payloads 

20  04  Aeroballistics 

19.  A8STRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

^This  report  describes  the  initial  phase  of  a  study  of  the  fluid  motion  in  a  con¬ 
tainer  that  is  spinning  and  nutating.  The  moments  exerted  by  this  fluid  motion 
on  the  container  have  been  demonstrated  to  cause  serious  flight  instabilities  of 
shells  with  liquid  payloads.  The  initial  phase  of  our  study  concentrates  on  the 
relevant  dimensionless  parameters  and  equations  that  govern  the  fluid  motion. 
Various  special  cases  are  discussed.  Our  analysis  and  aims  at  supporting 
ongoing  experiments  at  Chemical  Systems  Laboratory  (CSL).i*  and  Ballistics 
Research  Laboratory  (BRL)  and  at  providing  a  link  between  measurements 
of  despin  moments  and  side  moments.  It  also  prepares  for  an  in-depth  analysis 
of  the  interior  flow  field  for  comparison  with  flow  visualizations  at  CRDC. 

*Now  the  Chemical  Research  and  Development  Center  (CRDC),  US  Army  Armament, 
Munitions,  and  Chemical  Command. 


20  DISTRIBUTION /AVAILABILITY  OF  A8STRACT 
G3  'JNClASSIFiE D/U.NLIMITEO  □  SAME  AS  RPT 


22*  NAME  OF  RESPONSIBLE  INDIVIDUAL 

BREUDa  C.  ECKSTEIN 


DO  FORM  1473, 84MAR  83APRe 


ACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

as  rpt  G  otic  users  UNCLASSIFIED 


22b.  TELEPHONE  (Include  Area  Code)  22c.  OFnCE  SYMBOL 

(301)  671-2914  DRSMC-CLJ-IR  (A 


83  APR  edition  may  be  used  until  exhausted.  SECURITY  CLASSIFICATION  OF  this  PAGE 

Allother  editions  are  obsolete.  UNCLASSIFIED 


RlTY  CLASSl  FlC*  T.OH  0  F  THIS  PACE 


BLANK 


2 


SECURITY  Classification  C  *  Th;J  p*OE 


PREFACE 


The  work  described  in  this  report  was  authorized  under 
project  No.  1L161 102A7 1A ,  CB  Defense  Research.  This  work  was 
started  in  July  1982  and  completed  in  January  1983. 

The  use  of  trade  names  in  this  report  does  not  constitute 
an  official  endorsement  or  approval  of  the  use  of  such  commercial 
hardware  or  software.  This  report  may  not  be  cited  for  purposes  of 
advertisement . 

Reproduction  of  this  document  in  whole  or  in  part  is 
prohibited  except  with  permission  of  the  Commander,  Chemical 
Research  and  Development  Center,  ATTN:  DRSMC-CLJ-IR  (A),  Aberdeen 
Proving  Ground,  Maryland  21010-5423.  However,  the  Defense  Technical 
Information  Center  and  the  National  Technical  Information  Service 
are  authorized  to  reproduce  the  document  for  United  States 
government  purposes. 


Acknowledgments 


Major  parts  of  the  reported  work  were  supported  under  the 
Scientific  Services  Program,  Contract  No.  DAAG29-81-D-0100',  </'-: 

Delivery  Order  No.  0284  conducted  between  July  1982  and  January 
1983;  COTR  Miles  C.  Miller,  Chemical  Research  and  Development 
Center,  Aberdeen  Proving  Ground,  Maryland  21010-5423.  Additional 
work  was  supported  by  the  US  Army  Research  Office,  Mathematics 
Division,  under  Contract  DAAG29-82-K-01 29 . 


Accession  For 


,9 


9 


CONTENTS 


1.  INTRODUCTION  .  7 

2.  EXPERIMENTAL  APPROACH  AND  RESULTS  . .  9 

2.1  Despin  Moment  .  9 

2.2  Yaw  Angle  Growth  Rate  .  11 

2.3  Flow  Field  .  12 

3.  DIMENSIONAL  ANALYSIS  .  12 

3.1  Formal  Scaling  Aspects  .  18 

3.1.1  Scaled  Container  .  19 

3.1.2  Scaled  Velocities  o£  Working  Fluids  .  20 

3.1.3  A  Small  Model  Test  Fixture  .  20 

3.1.4  Data  Reduction  and  Cross-Checks  .  20 

3.2  Parametrical  Dependence  .  21 

3.2.1  An  Example:  Despin  Moments  22 .  22 

3.2.2  Gyroscope  Data  .  28 

4.  GOVERNING  EQUATIONS  .  29 

4.1  Coordinates  Fixed  to  the  Cylinder  .  29 

4.1.1  Dimensionless  Form  . 32 

4.1.2  Components  in  Cylindrical  Coordinates  .....  32 

4.2  Coordinates  Fixed  to  the  Axes  of  Rotation  .  34 

4.2.1  Decomposition  of  the  Flow  Field  .  35 

4.2.2  The  Reduced  Pressure  p'  .  37 

5.  LINEARIZED  EQUATIONS  . 38 

5.1  The  Case  of  Weak  Forcing  .  38 

5.2  The  Inviscid  Limit  .  39 

5.3  Effect  of  Viscosity  . 44 

5.4  Scaling  Aspects  .  45 

5.5  Estimates  on  the  Velocity  .  45 

6.  REMARKS  ON  MORE  GENERAL  CASES  . 47 


5 


LITERATURE  CITED 


51 


t 


0 


1. 

2. 

3. 


1. 


2. 


3. 


FIGURES 


Dimensional  Analysis  Test  Cylinder  Configuration  ....  13 

Deviation  from  the  Average  Value  Based  on 

Runs  74  and  25  .  24 

Modified  Plot  of  Data  from  Miller's 

Figure  1 36  27 

TABLES 

Despin  Moments  Analysis  Data — from  Figure  7  of 

Miller**  . 23 

2 

Values  of  M/fl  for  the  Four  6  Angles — 

from  Figures  8  and  9  of  Miller6  .  24 

Approximation  of  Miller's  Figure  12  with  the 

Cylinder  Filled  with  a  Homogeneous  Fluid  of 

High  Viscosity  .  25 


6 


FLUID  MOTION  IN  A  ROTATING  AND  NUTATING  CONTAINER 

1 .  INTRODUCTION 

Spin-stabilized  projectiles  with  liquid  payloads  can 

experience  severe  flight  instability.  Instrumented  flight  tests 

indicate  that  this  instability  is  characterized  by  a  rapid  increase 

in  projectile  yaw  angle  accompanied  by  an  abrupt  loss  in  spin  rate. 

Detailed  flight  data  for  the  XM761  smoke  screening  projectile  were 

12  3 

given  by  D'Amico.  '  '  These  data  clearly  show  that  the  instabil¬ 
ity  is  due  to  the  motion  of  liquid  white  phosphorus  embedded  in 
cotton  wicks;  at  lower  ambient  temperatures  when  the  phosphorus  is 
in  a  solid  state,  the  projectile  experiences  a  stable  flight. 

A  laboratory  test  fixture  was  developed  by  Miller  at  CSL* 

to  measure  the  despin  moment  of  a  full-scale  XM761  container  under- 

4  5  6 

going  spinning  and  nutation.  Miller's  tests  '  verified  the 
suggestion  of  Vaughn7  that  white  phosphorous  in  the  wicks  behaves 
like  a  homogeneous,  very  viscous  fluid.  The  tests  also  revealed  a 
maximum  despin  effect  for  fluids  of  kinematic  viscosity  in  the 
range  of  0.1  sq  m/sec.  Flight  tests  of  projectiles  having  a  con¬ 
tainer  filled  with  corn  syrup  of  kinematic  viscosity  0.2  sq  m/sec 
were  conducted  at  BRL  and  showed  instability  very  similar  to  that 
of  the  XM761. 

The  stability  problem  of  the  XM761  has,  meanwhile,  been 
overcome.  The  cotton  wicks  have  been  replaced  by  felt  wedges,  and 
separated  by  a  longitudinal  baffle  and  impermeable  plastic  foil. 

In  view  of  future  designs,  however,  there  is  ongoing  interest  in 
the  interior  fluid  motion  in  spinning  and  nutating  containers.  The 
experiments  at  CSL  have  significantly  simplified  the  research  tar¬ 
get  by  showing  that  the  basic  phenomena  can  be  studied  for  homo¬ 
genous  liquid  fills.  The  experiments  have  also  shown  that 
viscosity  plays  an  important,  if  not  dominating,  role  in  the 

*  Now  the  Chemical  Research  and  Development  Center,  US  Army 
Armament,  Munitions,  and  Chemical  Command. 
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problem  and  this  observation  makes  the  study  of  the  fluid  motion  a 
difficult  task.  Not  only  is  an  additional  parameter  introduced, 
but  moreover,  the  theoretical  analysis  of  the  motion  cannot  resort 
to  the  simpler  formulations  of  inviscid  theory  or  boundary-layer 
theory . 


Major  efforts  to  reveal  the  fluid  motion  and  to  determine 
the  despin  and  side  moments  required  for  aeroball istic  investiga¬ 
tions  have  been  purely  experimental.  Miller  conducted  void  obser- 

Q 

vations  in  a  transparent  cylinder.  He  also  developed  an  improved 
spin  fixture  that  allows  for  extended  despin  moment  measurements  as 
well  as  for  flow  visualizations  using  a  laser-induced  colored  dye 

Q 

technique.  D'Amico  and  Rogers  located  1 i qu id- f i 1 led  containers 
within  the  rotor  of  a  freely  gimballed  gyroscope  and  measured  yaw 
growth  rates  at  fixed  spin  rates. ^  The  yaw  growth  indicates  the 
side  moment  exerted  by  the  liquid  motion.  Observations  with  highly 
viscous  fluids  showed  that  the  growth  rates  were  inconsistent  with 
the  predictions  based  on  the  resonant  interaction  with  inertial 
oscillations  as  they  occur  in  low-viscosity  fluids.  '  Although 
the  measurements  of  despin  moments  and  yaw  growth  rates  lead  to 
consistent  conclusions  concerning  the  role  of  viscosity,  they  pro¬ 
vide  two  completely  separate  sets  of  data.  The  link  between  these 
sets  is  given  by  the  internal  flow  field  in  terms  of  velocity  and 
pressure  distribution.  This  flow  field,  however,  is  as  yet  unknown. 

A  first  step  into  linking  and  better  understanding  the 
available  data  is  to  identify  the  nondimensional  parameters  that 
model  the  flow.  The  experimental  data  base  is  as  yet  too  small  to 
definitely  identify  these  parameters,  in  particular  the  Reynolds 
number.  We,  therefore,  resort  to  formal  dimensional  analysis  and 
to  an  analysis,  of  governing  equations  and  special  solutions.  The 
availability  of  the  relevant  parameters  also  efficiently  supports 
the  experiments  since  the  amount  of  data  recording  can  be 
drastically  reduced  and  (future)  experiments  can  be  carried  out 
with  properly  scaled  containers,  spin  and  nutation  rates,  and 
scaled  working  fluids. 
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The  second  step  is  to  set  up  the  governing  equations  in 
the  best-suited  coordinate  system  and  to  prepare  for  approximate 
solutions.  Theoretical  work  in  this  direction  is  mostly  based  on 
the  Stewartson-Wedemeyer  approach,  using  an  inviscid  core  flow  that 
interacts  with  the  boundaries  through  viscous  boundary  layers. 
Essential  improvements  have  been  achieved  by  Murphy  and  Gerber, ^ 
Sedney  and  Bartos.  Use  of  the  boundary-layer  approximation, 
however,  precludes  application  of  these  theories  and  validity  of 
their  results  at  low  Reynolds  numbers.  The  restriction  to  small 
deviations  from  solid  body  rotation  is  always  necessary  for  working 
with  linearized  equations.  We  consider  various  sets  of  linearized 
equations  and  discuss  the  inherent  assumptions  in  relation  with 
estimates  from  experimental  data.  In  special  cases,  more  general 
solutions  or  new  insight  into  the  structure  of  the  motion  is 
obtained.  Finally,  we  discuss  some  aspects  of  computational 
methods  for  solving  the  complete  equations. 

2.  EXPERIMENTAL  APPROACH  AND  RESULTS 

Before  going  into  analysis,  we  briefly  describe  the 
experimental  methods,  the  range  of  parameters,  and  the  results 
obtained  to  date. 

2. 1  Despin  Moment. 

Measurements  of  the  despin  moment  were  conducted  using  a 

4 

test  fixture  at  CSL  described  by  Miller.  A  full-scale  XM761 
payload  canister  was  mounted  in  a  frame  under  an  angle  0  to  the 
vertical  axis.  The  canister  was  cylindrical,  with  inner  radius 
a  =  60.3  mm,  length  2c  =  517.6  mm,  and  a  total  mass  of  about 
12.25  kg.  The  angle  0  could  be  varied  to  be  5,  10,  15,  or 
2U  degrees.  The  cylinder  was  spun  up  by  an  air  turbine  to 
realistic  spin  rates  of  w  _<  6000  rpm.  Then,  the  frame  was  rotated 
aoout  its  vertical  axis  up  to  the  desired  nutation  rate  of 
fi  £  600  rpm.  At  higher  nutation  angles,  the  spin  rate  started  to 
decrease  before  the  frame  reached  high  nutation  rates  (Miller, 
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private  communication,  1982).  After  cutting  off  the  air  turbine, 
the  spin  rate  was  recorded  as  a  function  of  time.  The  total  despin 
moment  was  obtained  as  the  product  of  spin  deceleration  and  axial 
moment  of  inertia.  The  1  iqu id- induced  moment  was  found  by 
subtracting  the  :.ioment  due  to  friction  that  was  previously 
determined.  As  a  result,  the  1 iqu id- i nduced  despin  moment  is  given 
as  function  of  ft,  ^ ,  and  8.  The  method  implies  the  following 
assumptions : 

a.  The  nutation  rate  is  constant  during  spin-down. 

b.  The  motion  of  the  liquid  fill  has  no  effect  on  the 

moment  due  to  friction  as  a  function  of  ft,  to,  and  0. 

c.  The  fluid  motion  is  quasi-steady;  i.e.,  the  spin-down 
is  slow  enough  for  the  fluid  motion  to  fully  adapt  to 
the  instantaneous  spin  rates. 

d.  The  friction  moments  are  not  a  function  of  time. 

For  sufficiently  high  spin  rates  (a)  _>  1000  rpm),  the 
liquid  moment  was  found  to  be  independent  of  to.  Data  for  the 

XM761  payload  are  given  in  Miller's  figures  7  to  9.^  Additional 

data  were  obtained  for  the  cylindrical  container  filled  with 
homogeneous  liquids  of  kinematic  viscosities  in  the  range  from 
v  =  10  ^  Sq  m/sec  (water)  up  to  1.35  sq  m/sec  (corn  syrup) .  These 
data  are  shown  in  figures  12  to  14  of  the  same  paper.  No  attempt 
was  made  to  measure  the  side  moment  exerted  by  the  fluid,  which 
would  Lend  to  increase  the  nutation  (or  yaw)  angle  8  .  Only  one 
aspect  ratio,  c./a  =  4.29,  of  the  cylinder  was  studied  in  these 
experiments. 

An  improved  test  fixture  is  presently  in  operation  and 
more  detailed  data  will  soon  be  available. 


2.2 


Yaw  Angle  Growth  Rate . 


Measurements  of  the  yaw  angle  growth  rate  using  a 
gyroscope  whose  rotor  contained  a  liquid-filled  cylindrical  cavity 
were  conducted  by  Karpov, Scott  and  D'Amico*^,  and  D'Amico  and 
Rogers.*®  The  nutation  rate  of  the  (empty)  rotor  for  a  given  spin 
rate  was  varied  by  changing  the  moments  of  inertia.  The  rotor  was 
spun  up  with  the  axis  held  fixed,  and  time  was  allowed  for  the 
fluid  to  reach  a  state  of  solid-body  rotation.  After  releasing  the 
gyroscope,  the  spin  rate  was  held  constant.  The  yaw  angle  was 
measured  in  terms  of  an  amplitude,  A,  and  In  A  was  recorded  as  a 
function  of  time. 

The  earlier  experiments  in  the  range  of  high  Reynolds 
numbers  served  to  verify  Stewartson's  theory.  Linear  increase  of 
In  A  with  time  (constant  yaw  growth  rates)  was  found  in  some  range 
of  small  yaw  angles.  Scott  and  D'Amico  observed  nonlinearity  at 
angles  as  low  as  1  degree  (for  c/a  -  3).  Virtually,  a  changeover 
to  a  constant  growth  rate  occurred.  At  higher  viscosity,  the 
changeover  still  occurred  but  at  higher  yaw  angles.  The 
experiments  of  D'Amico  and  Rogers  were  made  with  silicone  oils  of 
high  viscosity  and  with  cylinders  of  aspect  ratios  c/a  =  3.126  and 
1.042.  Spin  rates  were  in  the  range  of  3000  rpm.  Constant  yaw 
angle  growth  rates  were  found  up  to  the  limit  of  the  apparatus  at 
5  degrees.  The  growth  rates  at  a  low  Reynolds  number  turned  out  to 
be  proportional  to  the  nutation  frequency  and,  therefore,  must  be 
due  to  a  mechanism  different  from  the  Stewartson-Wedemeyer  model. 
The  data  were  correlated  with  various  dimensionless  parameters; 
however,  a  simple  parametric  dependence  was  not  found. 

The  relation  of  the  yaw  angle  growth  rate  to  the  moments 
exerted  by  the  liquid  is  not  easily  obtainable.  The  moments  of 
inertia  are  not  given  for  the  various  data  points,  precluding 
determination  of  the  side  moment  via  the  coefficients  of  the 
stability  equation.  Neither  the  inplane  momen.  -he  despin 

moment  have  been  separately  measured.  There'.,  ct.  t-e  da  a  are 
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completely  separated  from  those  in  2.1.  They  lead,  however,  to  the 
same  conclusion,  that  a  basically  viscous  mechanism  can  lead  to  the 
instability  of  a  projectile  with  liquid  payload. 


2. 3  Flow  Field. 

The  flow  field  that  develops  from  the  viscous  fluid  motion 
is  unknown.  Miller  reported  void  observations  for  a  (95%)  filled 

O 

cylinder  with  an  aspect  ratio  of  4.45.  Due  to  the  high  spin  rate 
(u>  =  4000  rpm)  the  void  was  axisymmetric  for  solid-body  motion.  A 
characteristic  wavy  distortion  of  the  void  occurred  when: 

a.  The  nutation  rate  ft  was  increased  at  fixed  v  and  0, 

b.  The  nutation  angle  6  was  increased  at  fixed  v and  0, 

c.  The  kinematic  viscosity  was  decreased  at  fixed  0  and 

ft. 


For  low  viscosity,  the  void  was  irregular,  indicating  turbulent 
fluid  motion.  The  distortion  was  essentially  restricted  to  the 
plane  spanned  by  the  spin  axis  and  nutation  axis.  A  similar 
observation  was  sketched  by  Scott's  figure  1  without  giving  any 
detail.17  whereas  Scott  indicated  that  the  void  offset  from  the 
axis  at  the  cylinder  end  plates,  Miller's  photographs  show  no 
offset.  The  effect  of  the  void  on  the  fluid  motion  is  unknown. 

Flow  field  observations  are  presently  in  preparation  at 
CSL.  A  laser-induced  colored  dye  technique  combined  with  high¬ 
speed  photography  is  utilized  to  determine  the  three-dimensional 
velocity  field. 

3.  DIMENSIONAL  ANALYSIS 

We  consider  a  cylinder  of  length  2c  and  a  radius  a, 
rotating  about  its  axis  with  the  spin  rate  w.  The  axis  of  the 


cylinder  is  inclined  to  the  (vertical)  nutation  axis  by  the 
nutation  angle  6  and  rotates  with  nutation  rate  ft  about  the 
vertical.  The  cylinder  is  completely  filled  with  a  fluid  of 
constant  density  p  and  kinematic  viscosity  v.  The  cylinder  is 
symmetric  with  respect  to  the  margin  0  defined  by  the  point  of 
intersection  of  the  two  axes  of  rotation  (see  figure  1). 


spin  nutation 

axis  axis 

Figure  1.  Dimensional  Analysis  Test  Cylinder  Configuration 


The  motion  of  the  fluid  is  governed  by  the  continuity 


and  the  Navier-Stokes  equations 


DV  x  2 

-  =  -  i  yP  +  vr  v  (2) 

Dt  p  — 

where  V  is  the  velocity  vector  and  p  the  pressure.  The  V  in  equa¬ 
tion  (2)  is  measured  in  an  inertial  coordinate  system.  The  body 
force  due  to  gravity  g  is  neglected  in  equation  (2).  Justification 
will  be  given  shortly.  Equations  (1)  and  (2)  are  subject  to  the 
no-slip  and  no-penetration  conditions  at  the  side  wall  and  end 
walls  of  the  cylinder. 

The  solution  of  equations  (1)  and  (2)  under  the  given 
boundary  conditions  is  completely  determined  when  the  following 
quantities  are  given:  a,  c,  0,  uj,  fi,  p,  and  v.  Hence,  any  quantity 
derived  from  this  solution  can  be  written  as  a  function  of  these 
seven  parameters.  As  an  example,  we  write  the  despin  moment  M 
exerted  by  the  flow  at  the  cylinder  as 

M  =  t(a,  c ,  0,  w,  ft,  p,  X )  (3) 

The  principle  of  dimensional  homogeneity  requires  that  the  function 
of  f  has  the  dimension  of  the  moment  M.  According  to  Buckingham's 
Il-theorem,  this  can  only  be  satisfied  if  the  function  of  f  is 
composed  of  terms  that  are  products  of  powers  of  the  seven 
parameters.  °  For  example, 


f.  =  const  (a) 


(c) 


{  0) 


(col 


fi) 


( p)  ( v) 


P7 


(4) 


where  the  constant  is  dimensionless  and  Pj,...,P7  are  constant 
powers  to  be  determined.  We  utilize  the  following  dimension  table: 
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►. 


1 

I 


4 


Length  (L)  2  1  1  0  0  0  -3  2 

Time  (T)  -2  0  0  0  -1  -1  0  -1 

Mass  (m)  100  0  0  0  1  0 

2  2 

which  expresses,  for  example,  the  dimension  of  M  as  m  •  L  /T  . 
Exploiting  dimensional  homogeneity  of  equation  (4)  and  comparing 
powers  of  L,  T,  and  m  separately  provides  the  following  system  of 
equations  for  Pj,,,.,p7: 


L:  2  =  p1  +  p2 


(5a) 


a 

* 


! 


T: 


-2 


8  “P4  "  P5  " 


P  7 


(  5b) 


m:  1  =  p6  (5c) 

Since  only  three  equations  are  available  to  relate  the  seven 
quantities,  we  are  left  with  4  degrees  of  freedom  that  express: 

a.  The  presence  of  the  nondimensional  nutation  angle  6 

b.  The  presence  of  two  length  scales,  a  and  c 

c: .  The  presence  of  two  time  scales,  1/uj  and  1/ft 

d.  The  possibility  of  defining  an  additional  length 

scale  (for  example,  a2/v)  or  time  scale  (for  example, 
a/v/w v )  based  on  the  kinematic  viscosity. 

Only  equation  (5c)  provides  unambiguously  p^  =  1,  since  there  is  a 
unique  scale  for  the  mass.  Eliminating  p6  from  equation  (5a) 
results  in 


L:  5  =  Pl  +  p2  +2p? 


(5d) 


15 


In  order  to  remove  the  ambiguity,  we  choose  a,  uj,  and  p  to  form  a 
basic  system  of  dimensions.  We  can  then  express  equation  (4)  in 
the  form 


=  const 


p  +p2+2p7 

(  a )  1  1 


P4+P5+P7 

( til) 


where  the  powers  of  a ,  ui,  and  p  are  determined  by  equations  (5d), 
(5b),  and  (5c),  respectively.  Introducing  the  dimensionless  groups 

X  =  c/a  (aspect  ratio) 


t  *  ft/u>  (nutation  frequency) 
0  -  sin  0 
2 

Re  a  dia  /v  (Reynolds  number) 


(6) 


we  obtain 


2  5  /  P  2  P3 

f  ^  =  const  •  pur  a  yX  >0 

where  the  expression  in  parenthesis  is  dimensionless  and  the  powers 
P2'  P3'  P5'  andl  P7  are  arbitrary.  Forming  the  dimensionless  despin 
moment  M* ,  we  obtain  from  equations  (3)  and  (7)  in  general  form 

- j - r  =  M*  =  M*  (  Re  ,  X,  t,  0)  (8) 

puj  a3 

In  a  completely  analogous  manner,  we  can  form  other  dimensionless 
flow  quantities  with  the  functional  dependence  on  the  same 
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parameters  (see  equation  6).  Different  ways  to  find  the  functional 
dependence  will  be  discussed  in  Section  3.2. 

If  partially  filled  cylinders  were  to  be  considered,  the 

list  of  parameters  (equation  6)  had  to  be  completed  by  the  fill 

ratio  =  (fluid  volume )/( volume  of  the  cylinder).  If  gravity,  g, 

2 

were  taken  into  account,  the  Froude  number,  Fr  =  u>  a/g,  would  have 
to  be  included  in  the  list  of  parameters.  In  our  application,  Fr 
is  typically  of  the  ordei  10  — the  inertial  and/or  viscous  forces 
dominate  the  fluid  motion.  The  negligible  effect  of  the  Froude 
number  is  clearly  demonstrated  by  the  cylindrical  free  surface  in 

g 

Miller's  void  observations  (see  figures  3a,  4a,  and  5d) .  This 
surface  should  be  an  axisymmetric  paraboloid,  if  gravity  had  an 
effect.  Anyway,  if  the  fluid  boundaries  are  fixed,  as  in  the  case 
of  a  completely  filled  cylinder,  the  Froude  number  is  insignifi¬ 
cant,  and  the  only  effect  of  gravity  is  to  contribute  a  hydrostatic 
distribution  to  the  pressure. 

The  definition  of  the  Reynolds  number  in  equation  (6)  is 
as  arbitrary  as  the  choice  of  the  basic  scales  a  and  u>.  Other 
Reynolds  numbers  can  be  derived: 

Re  ^  =  Re  •  T  ,  Re  2  =  Re  •  t  •  1 2  r  Re  3  =  Re  •  T  (  X  o )  2  (9) 

or  in  terms  of  physical  parameters 

Re ^  =  fia2/v,  Re ^  =fic2/v,  Re3  =fib2/v  (10) 

where  b  =  c  •  sin  8  is  the  nutation  radius  of  the  cylinder  (see 
figure  1).  The  question,  which  of  these  Reynolds  numbers  is  most 
convenient,  that  is,  leads  to  the  simplest  functional  relation  for 
M*  and  so  on,  cannot  be  answered  by  the  present  forma1  dimensional 
analys  is . 
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3.1 


Formal  Scaling  Aspects 


There  are  two  important  consequences  of  this  dimensional 
analysis,  however.  First,  we  have  reduced  the  number  of  seven 
physical  parameters  by  the  number  of  three  basic  dimensions  to  only 
four  dimensionless  parameters.  This  simplifies  the  mapping  of  the 
parametrical  dependence,  second,  equations  such  as  (8)  show  that 
the  nondimensional  flow  quantities  are  identical  if  Re,  X,  x,  and  a 
are  identical. ,  no  matter  how  these  values  are  achieved.  This  fact 
allows  working  with  properly  scaled  physical  quantities. 


Let  us  replace  the  original  parameters  a,  c,...,  by  the 
model  parameters  am  =  a  •  ota ,  cm  =  c  *  <*c,...,  w^ere  <*=,»  ac '  '  '  *  ' 
are  scale  factors.  Then,  a  physically  similar  solution  (identical 
in  nondimensional  terms)  is  obtained  if  the  four  parameters 
(equation  6)  remain  unchanged.  In  other  words,  if 


a  =  a  ( l°a> 

C  3 

a  =  a  (10b) 

ft  co 

a  •  =  a  (10c) 

co  a  v 

a-  =  1  (10a) 

y 

This  last  equation  shows  that  the  nutation  angle  (as  £ 
dimensionless  quantity)  cannot  be  scaled.  Witn  equation  (10) 
satisfied,  the  same  value  of  M*  =  M*  will  be  obtained  according  to 
equation  (8).  The  physical  (dimensional)  moment,  however,  will 
change  as 


M  = 
m 


a  pa 


co  a  a 


M 


(  11  ) 


The  relations  (equutions  10  and  11)  allow  model  tests  without  any 
more  detailed  knowledge  of  the  physics  of  the  problem.  It  is 
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mainly  equation  (10c)  that  allows  for  changes  of  rotation  rates, 
size  of  the  cylinder,  and  viscosity  of  the  liquid. 

3.1.1  Scaled  Container. 

The  experiments  reviewed  in  Section  2.1  typically  use  a 
full-scale  container  with  a  mass  of  12.25  kg,  rotate  it  at  5000  rpm 
(10,000  rpm  in  the  new  test  fixture),  and  impose  a  nutation  with 
500  rpm  under  an  angle  of  up  to  20  degrees.  Clearly,  operating 
this  experiment  bears  the  danger  of  catastrophic  failure  and, 
therefore,  requires  extreme  care  in  designing,  machining,  and 
balancing  the  apparatus.  The  reason  is  obvious  from  equation  (7) 
or  (11):  the  moments  increase  with  the  fifth  power  of  the  length 
and  with  the  second  power  of  the  rotation  rates.  The  scaling 
relations  (equation  10)  allow  the  design  of  a  harmless,  smaller 
test  fixture  at  reduced  cost  that  bears  other  advantages. 

Let  us  first  consider  that  the  working  fluid  remains 

unchanged,  (ap  =  1,  =  1 ) .  If  the  linear  dimensions  of  the 

cylinder  are  reduced  by  a  factor  of  a  <  1,  equation  (10c) 

2  ** 

immediately  requires  a,  =  1/a  ,  an  increase  of  the  rotation  rates 
9  co  a 

by  a  factor  1/or.  From  equation  (11)  it  becomes  clear  that  the 
a 

moments  are  only  reduced  by  a  factor  aa .  A  drastic  reduction  of 

the  moments  cannot  be  achieved  in  this  manner,  since  the  required 

rotation  rates  would  be  prohibitive.  It  seems  necessary, 

therefore,  to  change  the  working  fluid.  For  liquids,  there  will  be 

only  a  small  variation  in  density,  ap  =  1,  and  the  variation  will 

be  mainly  in  the  viscosity.  Keeping  the  original  rotation  rates, 

2 

a  *  1  (equation  10c)  requires  a  =  a  Using  a  model  with  «  = 

(x)  V  d  a 

1/2,  say,  a  liquid  with  vm  =  v/4  must  be  used  to  maintain  dynamical 
similarity.  In  view  of  the  high  viscosities  of  interest, 
such  a  liquid  can  easily  be  found.  The  moments,  then,  are 
considerably  reduced  by  a  factor  a  a  -  1/32.  Further  reductions 

P  a 

can  be  achieved  by  the  use  of  even  smaller  models. 
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3.1.2 


Scaled  Velocities  o f  Working  Fluids. 


Disregarding  structural  problems,  scaling  may  have  other 

benefits,  (for  example,  flow  visualizations).  Should  it  turn  out 

that  the  velocities  are  too  large  for  easy  visualization  of  the 

colored  dye,  the  time  scale  of  the  phenomena  in  the  full-size 

cylinder  can  be  changed  by  reducing  spin  rate,  nutation  rate,  and 

viscosity  by  the  same  factor,  =  a^.  The  moments  drop  by 

a  a2 ,  slowing  down  the  spin-down  process.  The  velocities  drop  by 

p  oo  3 

the  factor  a  .  Smaller  picture  frequency  (factor  ay)  and  larger 
exposure  times  (factor  1/a  )  would  recover  the  original  results. 
Should  it  be  desirable  to  work  with  a  special  fluid,  for  example, 
water  instead  of  some  more  viscous  silicone  oil,  the  necessary 
adjustments  of  rotation  rates  and  cylinder  size  can  be  made 
according  to  a,  =  a  a..,  from  equation  (10c). 

V  o  (JU 


3.1.3 


A  Small  Model  Test  Fixture. 


By  proper  choice  of  the  scale!  factor  a,  it  seems  possible 
to  set  up  a  model  test  fixture  (based  on  a  standard  record  player 
as  a  turntable)  at  low  cost  for  direct  visual  observation  of  the 
motion  in  high-viscosity  fluids.  With  fixed  =  45  rpm,  the  range 
of  nutation  rates  can  be  covered  with  a ^  >  0.075.  A  reasonable 
range  of  x  can  be  obtained  with  o>m  <  2000.  With  aa  =  ac  “  0.4,  the 
filled  model  cylinder  has  a  mass  of  approximately  0.8  kg.  Moments 
are  reduced  by  about  tour  orders  of  magnitude.  The  resulting 
factors  a  >  0.012  allow  covering  the  interesting  range  of  higher 
viscosities. 

3.1.4  Data  Reduction  and  Cross-Checks. 


Miller's  experiments  at  CSL  provide  despin  moments  for 
fixed  X  =  c/a,  for  four  values  of  o  =  sin  0,  and  numerous  values  of 
the  viscosity.  Nutation  rates  can  be  varied  in  a  broad  range,  and 
spin  rates  vary  continuously  during  spin-down.  It  is  obvious  from 
equation  (8)  that  only  two  parameters,  Re  and  x,  need  to  be  varied 
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to  obtain  (or  document)  the  complete  parametric  dependence  of  M* 
for  fixed  X  and  o.  It  is  useful,  therefore,  to  present  all  data  in 
dimensionless  form.  The  experimental  procedure  mentioned  earlier 
partially  duplicates  data.  As  an  example,  nutation  rates  and 
^2  =  with  fluids  of  viscosities  and  \>2  =  v^/2  should 

provide  identical  values  of  M*  at  end  =  wi/2.  This 
duplication  can  be  utilized  for  a  cross-check  of  the  data  and  for 
verifying  the  validity  of  the  assumptions  outlined  in  Section  2,1. 
With  this  validity  assured,  the  duplication  can  be  avoided  in 
future  experiments  by  using  cylinders  of  different  aspect  ratios. 

3 . 2  Parametrical  Dependence. 

The  second  step  in  the  fluid  motion  analysis  and  the 
resulting  moments  requires  finding  the  functional  dependencies  of 
nondimensional  variables  such  as  M*  in  equation  (8)  on  the  param¬ 
eters  Re,  X,  T,  and  a.  The  results  can  then  be  recast  in  terms  of 
other,  perhaps  better  suited  parameters,  such  as  one  of  the  modi¬ 
fied  Reynolds  numbers  in  equations  (9)  and  (10).  Various  methods 
can  be  exploited  and  their  results  combined  to  achieve  this  goal: 

m  physical  ai-yuuicntS 

•  Analysis  of  experimental  data 

•  Formal  analysis  of  the  governing  equations 

•  Analysis  of  numerical  results 

•  Analytical  solutions. 

Knowledge  of  the  analytical  solution  for  the  velocity  and  pressure 
field  would  completely  solve  the  problem.  However,  there  is  no 
chance  to  tind  the  solution  of  the  full  Navier-Stokes  equations. 
Simplifications  that  may  permit  findiny  approximate  solutions  will 
be  discussed  in  Section  5.  Numerical  solutions  of  the  full 
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equations  can  be  obtained,  in  principle,  but  are  not  yet  available. 
Analysis  of  the  governing  equat ions — even  without  solving  them — 
sometimes  reveals  special  combinations  of  the  parameters  as  most 
relevant.  An  example  is  the  Taylor  number  Ta  =  2Re  .  ( -  R^)/ 

( R2  +  Rj),  with  Re  =  f2  R  ^  ( R  2  ~  R  ^)/v,  that  characterizes  the  fluid 
motion  of  a  viscous  liquid  between  concentric  cylinders  of  radii 
and  R  ^  ,  with  the  outer  fixed  and  the  inner  rotating  with  ft.  Simi¬ 
lar  combinations,  square  of  the  Reynolds  number  times  a  purely 
geometric  parameter,  occur  in  other  problems.  Physical  arguments 
are  sparse  in  the  present  context  and  can  be  misleading,  in  fact, 
since  the  motions  of  rotating  fluids  are  often  beyond  our 
imagination.  In  the  following  discussion,  we  attempt  an  analysis 
of  the  experimental  data, 

3.2.1  An  Example:  Despin  Moments. 

This  analysis  is  heavily  restricted  by  the  small  number  of 
available  data.  These  data  are  uncertain  due  to  experimental 
errors  that  are  yet  known.  (Estimates'  from  the  cross-checks 
suggested  in  Subsection  3.1.4  would  be  helpful.)  Additional  error 
is  introduced  by  reading  the  data  from  small  figures. 

From  Miller's  figure  7,  it  appears  that  the  despin  moment 

2  6 

M  at  fixed  X  and  0  varies  proportionally  to  ft  .  Liquid  density 

and  viscosity  are  not  easy  to  identify  for  the  wick-type  fill.  We 

have  therefore,  not  used  M*  but  calculated  M/ft^.  The  values  of  t 
2 

and  of  M/ft  are  given  in  table  1,  There  is  only  a  small  scatter  in 

2 

the  values  of  M/ft  ,  largest  for  the  data  points  with  small  values 
of  M  that  are  difficult  to  read  accurately  from  the  graphs.  An 
average  value  of  1.77  x  10  Nms  is  obtained  such  that 

M  =  1.77  x  10~3  ft2  or  V-  =  1.77  x  10~3  T2  (12) 

Ou 

where  ft  or  00  are  in  radians/sec.  Based  on  runs  74  ana  2b  only, 

-  3  2 

the  average  value  would  be  smaller,  1.71  x  10  N  •  m  •  sec  . 
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Table 


Despin  Moments  Analysis  Data 


★ 


Run 

0) 

ft 

T 

M  10 J  M/ft2 

( rpm) 

(  rpm 

t 

<N  •  m)  (N  • 

2  , 

m  •  sec  ) 

74 

1000 

381 

0.381 

2.56 

1.61 

1250 

379 

0.303 

2 . 69 

1 .71 

1500 

366 

0.244 

2,70 

1.84 

2000 

347 

0. 174 

2.33 

1.76 

3000 

299 

0.0997 

1.60 

1.63 

4000 

225 

0.0564 

0.895 

1.61 

4500 

159 

0.0354 

0.479 

1.73 

25 

1500 

354 

0.236 

2.34 

1 .  70 

2000 

339 

0. 170 

2.26 

1.79 

3000 

293 

0.0977 

1.61 

1 .71 

4000 

221 

0.0553 

0.947 

1.77 

23 

2500 

296 

0.118 

1  .71 

1 . 78 

4000 

192 

0.0479 

0.652 

1.61 

4500 

100 

0.0222 

0.289 

2.64 

26 

2000 

192 

0.0958 

0.698 

1.73 

2500 

107 

0.0427 

0.277 

2.21 

21 

2000 

156 

0.0780 

0.393 

1.47 

2500 

101 

0.0403 

0.173 

1 . 55 

•Extracted  from 

Miller' 

s  Figure  7. 6 

The  deviation  from  this 

average  is 

shown  in  figure  2. 

It  is  diffi- 

cult  to  decide 

on  a  systematic  variation  with  u>  [or  the  Reynolds 

number  ( Re ) ] . 

Although 

there  seems 

to  be  a  tendency 

for  the  data 

of  these 

two  runs,  we  refrain  from 

any  additional  conclusion. 

From  Killer's 

figures  8  and  9,  the  effect  of 

6  or  o  can  be 

estimated 

.  The 

values 

of  M/ft2  for 

the  four  angles  6 

arc  given  in 

table  2. 

It  turns  out 

that  the  ratio  of  M./(ft  sin  6)  2 

is  almost  the 

same  for 

each  of  the  four  angles  . 

Significant  deviation  occurs 

only  for 

6  =  5 

deg rees  . 

For  6  >  5 

degrees,  the  data 

can  be 

approximated  by 


M  -  0 . 0 1 56  ( ft  sin  8)2  or  =  0.0156(to)2 

u> 


(  12) 
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a  run  74 
o  run  25 


Figure  2.  Deviation  from  the  Average  Value  Based  on  Runs  74  and  25 
Table  2.  Values  o£  M/ft2  for  the  Four  9  Angles* 


=  s in  0 


103  M/ft2 


M/  (ft  sin  0) 


Average 


(  degrees ) 

(N  •  m  •  sec  ) 

( N  •  m  ■  sec  ) 

5 

0.0872 

0.076 

0.0100  j 

10 

0.1736 

0,497 

0.0165  ( 

)  0.0156 

15 

0.2588 

1.02 

0.0153  ( 

20 

0.3420 

1  .  77 

0.151  ! 

From  Figures  8  and  9  of  Miller 


A  similar  approximation  can  be  obtained  from  the  data  in  table  3 
for  Miller's  figure  12.  In  this  case,  the  cylinder  is  filled  with 
a  homogeneous  fluid  of  high  viscosity.  Combining  data  for  all  four 
angles  0  results  in 

M  =  0.0u6i4(fi  sin  6)2  (14) 

The  different  constants  in  equations  (13)  and  (14)  can  be  attrib¬ 
uted  to  various  effects:  different  mass  of  liquid,  different 
effective  viscosity,  and,  in  addition,  different  geometry  due  to 
the  baffle  in  the  wick-type  payload. 

Table  3.  Approximation  of  Miller's6  Data  with  the  Cylinder 
Filled  with  a  Homogeneous  Fluid  of  High  Viscosity* 


0 

o  =  sin  0 

103  M/Q 2 

M/(Q  sin  0)2 

Average 

( degrees ) 

( N  •  m  •  sec 2 ) 

2 

( N  •  m  •  sec  ) 

5 

0.0872 

0.0689 

0.00904  \ 

10 

0.1736 

0.2116 

0.0070  ! 

0  00814 

15 

0.2588 

0.5559 

0.00830  ( 

20 

0.3420 

0.9581 

0.00819  ) 

•Miller's  Figure  126 


The  expressions  for  M  contain  only  quantities  that  are 
related  to  nutation.  This  reflects  Miller's  observation  that  the 
despin  moment  was  essentially  independent  of  the  spin  rate  when¬ 
ever  '■>  was  sufficiently  large.  The  other  parameters,  in  particular 
those  that  were  held  constant,  are  hidden  in  the  constants.  The 
aim  of  finding  the  functional  relation  [see  equation  (8)]  in  an 
as-simple-as-possiole  form  also  requires  reviewing  the  reference 
quantities  that  were  chosen  to  obtain  the  parameters  (see  equa¬ 
tion  (6)).  tquauions  (13)  and  (14)  suggest  that  the  reference 
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velocity  ma  should  be  replaced  by  ftc  sin  6.  From  physical  arcu- 

3 

ments ,  it  also  seems  reasonable  to  replace  pa  by  the  effective 

•5 

mass  of  liquid,  =  2Trpa“c,  so  that 

M  =  m£  (ftc  sin  6)2  •  M*  (15) 

3  2  2 

where  =  M*/2(irX  T  a  )  is  almost  independent  of  ft,  w,  and  0. 

The  dependence  of  on  the  viscosity  can  be  estimated 

from  Miller's  figure  13.^  Figure  3  shows  a  modified  plot  of  these 

data  in  the  form,  log  M  versus  log  1/v.  This  representation  has 

been  chosen  since  v  appears  only  in  the  Reynolds  number  and, 

therefore,  the  dependence  of  M*  on  Re  should  be  revealed  by  this 

plot.  Although  only  a  few  data  points  are  available,  it  is  obvious 

that  at  least  two  regions  can  be  distinguished.  For  high  viscosity 

2 

(region  I),  v  _>  0.2  m  /sec,  it  seems  that  M^  ~  1/v,  whereas  for 
v  <  0.1  m2/sec  (region  II),  the  data  align  along  Mj  -  ( 1 / A )  ®*^2, 
with  some  transition  region  in  between.  It  is  not  clear,  though, 
whether  a  second  transition  occurs  for  v  <  10  mVsec. 

The  data  of  Miller's  figure  12  and  the  moments  according 
to  equation  (14)  are  located  in  the  region  where  M|  ~  1/v.  Since 
ft  and  6  are  constant  for  figure  13  and  v  is  constant  for  figure  12, 
both  sets  of  data  can  be  combined  into 

2 

M  -  m.ft—  -  10  —  M-  =  m  0  ft  Re  ,  M  (16) 

X,  V  t  %  i  t 

where  M2  is  almost  independent  of  ft,  u>,  6,  and  v  From  comparison 
with  equation  (14)  we  obtain 

M2  =  0.00294  m2/sec  (17) 

The  value  of  the  Reynolds  number  for  the  data  in  figure  12  is 
Re^  -  2.05,  very  small  indeed.  For  the  range  of  smaller 
viscosities,  expressions  analogous  to  equation  (16)  cannot  be 
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derived  due  to  the  lack  of  data.  We  also  note  that  the 
approximations  suggested  in  this  section  need  further  verification. 

It  one  adopts  Re  3  as  the  relevant  Reynolds  number,  the 

void  observations  of  Miller  can  be  linked  to  the  despin  moments  in 

8 

figure  3.  The  slight  difference  in  the  lengths  of  the  cylinders 
can  be  neglected.  The  arrows  in  figure  3  indicate  the  Reynolds 
numbers  for  the  photographs  in  Miller's  figures  3,  4,  and  5.  To 
within  the  accuracy  of  reading  the  void  distortion  from  the  small 
reproductions,  the  value  of  Re3  sorts  the  pictures  in  the  sequence 
of  increasing  amplitude  of  the  sinusoidal  distortion.  The  void 
distortion  is  barely  visible  in  pictures  4a  and  5b,  but  increases 
as  Re3  changes  through  the  transition  region  into  region  II.  It  is 
tempting,  therefore,  to  associate  the  void  distortion  with  the 
deviation  of  the  despin  moments  from  equation  (16)  for  region  I. 

One  might  also  speculate  that  the  motion  in  region  I  is  essentially 
axisymmetr ic .  Flow  visualizations  of  the  axial  velocity  component 
can  decide  this  issue. 


3.2.2 


Gyroscope  Data. 


We  attempted  to  exploit  the  measurements  of  D'Amico  and 

Rogers  in  order  to  check  the  previous  conclusions  and  to  obtain 

additional  insight.10  However,  this  attempt  failed  since  the 

inertial  moments  of  the  gyroscope  are  not  given  in  the  paper;  these 

data  are  necessary  for  extracting  the  liquid  side  moments  from  the 

1  9 

equations  provided  by  the  tri-cyclic  theory. 


We  also  attempted  to  correlate  the  yaw  angle  growth  rates 
with  various  dimensionless  parameters.  One  of  our  observations 
seems  worth  reporting:  If  ^/Re  is  plotted  versus  Tr,  the  (few) 
data  points  fit  well  on  straight  lines,  in  particular  for  the 
cylinder  with  X  =1.  The  slope  of  these  lines  seems  to  be 
independent  of  w,  but  increases  with  viscosity. 


4. 


GOVERNING  EQUATIONS 


Since  the  experimental  data  base  is  too  small  tor  a 
detailed  analysis  of  the  moments  exerted  by  the  moving  fluid,  we 
prepare  for  complementary  theoretical  studies.  This  requires  one 
to  decide  on  a  suitable  coordinate  system  and  to  write  the 
equations  of  fluid  motion  and  appropriate  boundary  conditions  in 
these  coordinates.  Subsequently,  simplified  sets  of  equations  can 
be  developed  using  various  assumptions  for  approximate  solutions. 

The  geometry  and  kinematics  of  our  problem  suggest  one  of 
the  following  three  coordinate  systems: 

•  Inertial  or  earth-fixed  axes,  index  i 

•  Cyl inder- fixed  or  body-fixed  axes,  index  b 

•  Intermediate,  nutation-oriented  coordinate  system, 
index  n. 

The  choice  of  a  special  system  is  largely  guided  by  the  idea  that 
the  parameters  of  the  problem  should  appear  in  the  differential 
equations,  not  in  the  boundary  conditions.  The  boundary  conditions 
should  be  as  simple  as  possible.  These  requirements  are  obviously 
not  satisfied  if  equations  (1)  and  (2)  for  an  inertial  system  are 
used.  Only  v  and  p  appear  in  these  equations,  while  the  remaining 
five  parameters  are  hidden  in  the  boundary  conditions.  Moreover, 
the  boundary  conditions  are  time-dependent. 

4  .  1  Coordinates  Fixed  to  the  Cylinder . 

Coordinates  fixed  to  the  spinning  cylinder  are  analogous 
to  the  body  axes  often  used  in  aeroba 1 1  i st ics .  The  velocity  V  = 
in  equations  (1)  and  (2)  is  muasured  with  respect  to  an  inertial 
system.  The  body-fixed  system  rotates  with  angular  velocity  n  , 

where 


t )  =  jq  +  u(t) 


(  18) 


with  respect  to  the  inertial  system,  where  n  and  oj  point  in  the 
direction  of  the  nutation  axis  and  spin  axis,  respectively.  In 
this  system,  a  different  velocity,  will  be  measured. 

Velocities  yi  and  Vb  are  related  by 

Vi  -  vb  +  Qb  x  r 
20 

where  r  is  the  position  vector.  This  transformation 
continuity  equation  unchanged: 

v  •  Vj,  «  o  (20) 

Moreover,  for  a  Newtonian  fluid,  the  viscous  terms  are  invariant 
under  the  transformation  equation  (19).  Major  changes  occur  only 
in  the  acceleration  terms: 


(19) 

leaves  the 


DV  DYb 

dT  '  5T  +  %  *  +  *  ‘fib  x  +  V  E  <21) 

where  flb  denotes  the  time-derivative  of  2b  •  The  three  additional 
terms  containing  fib  are  the  Coriolis  acceleration,  the  centripetal 
acceleration,  and  the  acceleration  due  to  the  change  of  the 
rotation  rate,  respectively.  The  centripetal  acceleration  can  be 
written  in  terms  of  a  potential  function  <$c; 

£b  x  %  x  I>  “  -V*c  (22) 


and  can  be  considered  as  an  additional  (conservative)  force  per 
unit  :ic  3r  ■  Therefore,  this  term  is  written  on  the  right-hand  side 
of  the  momentum  equation: 


DVb  1 

5T  +  2nb  x  Vb  +  Pb  x  r  =  -  -p  Vp  + 


V*  +  nv  v. 

c  — b 


(23) 
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It  is  obvious  from  equation  (23)  that  the  centrifugal  force  can  be 
combined  with  the  pressure  p  into  a  reduced  pressure: 

Pb  =  P  *  P*c  =  P  -  \%  x  r  | 2  (24) 

The  boundary  conditions  of  no-penetration  and  no-slip  reduce  to  the 
simple  form 


Vb  s  o  at  the  boundaries  (25) 

Only  the  radius  a  and  the  half-length  c  of  the  cylinder  are  intro¬ 
duced  by  equation  (25).  All  the  other  parameters  of  the  problem 
are  contained  in  equation  (23)  with  ftb  from  equation  (18).  The 
system  of  equations  (20)  and  (23),  with  (25),  supports  trivial 
solutions  Vb  -  0,  pb  =  0,  whenever  Qb  =  0  (that  is,  rigid-body 
motion  of  the  fluid).  In  a  more  obvious  form,  this  result  can  be 
written  as 

Yh  =  0,  pb  i  0  for  wfi  sin  0*0  (26) 

There  are,  in  fact,  three  separate  cases  of  rigid-body  motion: 

•  For  zero  spin  rate,  w  =  0 

•  For  zero  nutation  rate,  P  =  0 

•  For  a  rotation  with  o  +  n  about  the  same  axis. 

According  to  equation  (24)  pb  =  0  describes  the  pressure 
distribution 


where  rd  is  the  distance  from  the  axis  of  rotation. 


(  27 ) 
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A  fourth  case  of  rigid-body  motion  occurs  for  large  vis- 

2 

cosity,  u  ■+•  06.  in  this  limit,  equation  (23)  reduces  to  v  Vb  =  0 
for  Stokes  flow. 

4.1,1  Dimensionless  Form* 

Using  the  basic  system  of  dimensions  (a,  w,  and  p),  we 
introduce  dimensionless  variables  (denoted  by  * )  in  the  following 
way : 


£  =  r  •  a ,  t  =  t/u) 

V  =  V  •  aw,  =  Pb  =  Pb  ‘  P<a<*>)2 

The  basic  equations  then  take  the  following  forms 


V 


*  ^b 


£  =  -7p. 


V2  0 


Re 


-b 


(28) 

(29) 


(30) 

(31) 


In  equation  (31),  7  refers  to  dimensionless  variables. 

We  note  that  the  spin  rate  co  must  be  nonzero  for  the 
availability  as  a  reference  quantity.  This  is  no  restriction,  how¬ 
ever,  since  according  to  equation  (26)  rigid-body  motion  occurs  for 
w  =  0.  We  are  rather  interested  in  the  deviations  from  rigid-body 
motions  for  w  /  0. 

4.1.2  Components  in  Cylindrical  Coordinates. 

In  view  of  the  axisymmetry  of  the  boundary,  we  introduce 
cylinder  coordinates  r,  4> ,  and  z,  where  the  z-axis  coincides  with 
the  cylinder  axis  and  r=l,  <f>  =  z  =  0  describes  a  fixed  point  of 
the  cylindrical  side  wall.  Since  no  further  use  is  made  of 
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where 


D' 


and 

D" 


3_ 

3r 


r  3? 


u 


z 


3  z 


2  2  2 

3  ,13  .1  3^  ,  d 

T^2  r  3r  r2  3$2  3z2 


The  boundary  conditions  require 


u  =  u.  =  u  =  0  for  r  =  1,  |  z  |  £  X,  0  £  $  <  2  tt  (35a) 

f  T  Z 

ur  =  u^  =  uz  =  0  for  z  =  +  X,  0  £  r  £  1,  0  <_  <p  <  2ir  (35b) 

Although  these  conditions  are  independent  of  time,  a  time-dependent 
solution  is  desirable,  due  to  the  time-dependence  of  the  coeffi¬ 
cients  n r (see  equation  (32)]  .  It  is  obvious  that  a  simple 
time-periodicity  with  period  T  =  2ir  will  evolve  if  transient  behav¬ 
ior  (after  sudden  start  of  nutation,  for  example)  is  disregarded. 


The  possibility  of  studying  small  deviations  from  rigid- 
body  motion  as  a  perturbation  of  the  zero-state  [see  equation  (26)] 
seems  to  support  the  use  of  this  body-fixed  coordinate  system.  The 
periodicity  in  time  can  be  taken  into  account  by  introducing  a 
modified  azimuthal  coordinate,  4> '  =  $  +  t .  At  closer  analysis, 
however,  this  change  of  coordinates  is  equivalent  to  working  in  a 
system  that  rotates  only  with  fi  about  the  nutation  axis. 

4.2  Coordinates  Fixed  to  the  Axes  of  Rotation. 


This  nutating  coordinate  system  rotates  with  constant 
angular  velocity  Qn  =  about  the  origin.  The  equations  of 
Section  4.1  require  little  change:  Quantities  with  index  b  are 
replaced  by  quantities  with  index  n  and  the  terms  with  fin  are 
dropped.  The  boundary  conditions,  however,  are  inhomogeneous,  due 
to  the  rotation  rate  w  of  the  cylinder  with  respect  to  the  nutating 
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system.  Instead  of  equation  (26),  the  trivial  solutions  for 
ft  sin  0=0  are  rigid-body  rotations  with  tu  about  the  cylinder 
axis.  Homogeneous  boundary  conditions  can  be  retrieved,  and  an 
apparently  simple  formulation  can  be  obtained  by  splitting  the 
velocity  field  into  a  rigid-body  rotation  with  w  and  a  deviation 
due  to  Q  sin  0/0.  In  the  following,  we  return  to  dimensionless 
variables . 


We  introduce  cylindrical  coordinates  r,  <p ,  and  z,  where 
the  z-axis  coincides  with  the  cylinder  axis,  as  before.  However, 
=  0  describes  points  in  the  plane  spanned  by  the  nutation  axis 
and  spin  axis,  because  0/0.  The  components  of  Q  =  ft  are 


ft  r  3  -  to  cos  #  , 


ft^  =  to  s i n  4>  ,  pz  =  t  ( 1  -  o2)1//2 


(36) 


Hence,  in  the  nutating  frame,  the  Coriolis  acceleration  introduces 
no  explicit  time-dependence. 

4.2.1  Decomposition  of  the  Flow  Field. 

We  decompose  the  velocity  Vn  =  (ur,  ,  uz  )  and  the 
reduced  pressure  pn  according  to 


u„  =  V.  UA  =  ft.r  +  v  u 
r  r  r  (£  <J> '  •  z  z 


V.,  U_  =  V. 


Pn  ~  P '  +  \  (1  +  2ftz)  r2 


(37) 


Although  the  continuity  equation  remains  in  the  form 


1  JS_ 
v  9r 


(rv  ) 
r 


+ 


1 

r 


0 


(  38) 


for  the  deviation  (v  ,  v,  ,  v  ),  the  j„ omentum  equations  take  the 

IT  y  Z 

following  form: 
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(39a) 


D'vr  -  4  -  2(1  +  n2>  v*  .  2«<,v2 


D'v*  +  “r1  *  2(1  *  nz]  vr  '  2nrvz 


i  l£l  +  J_  [D»„  .  \  .  L.  21e] 

r  3$  Re  [  <p  r2  r2  dp  J 


D'v  +  2ft  v, 
z  c  p 


P  r 


-  -r  -  2^r  *  Is  »■*. 


D'  now  reads: 


(39b) 


(39c) 


D' 


3 

Jc 


±_ 

dp 


Vz  Tz 


The  boundary  conditions  are  homogeneous: 


(40) 


Vr  =  v4»  =  vz  "  0 

for  r  *  1,  |z|  <  x  ,  0  <  <p  <  2n  (41a) 

vr  =  v4>  “  vz  “  0 

for  z  =  +  A ,  0£r<_l,0<_(f><2iT  (41b) 

The  system  equations  (38)  through  (41)  for  the  deviation  from 
rigid-body  motion  bears  some  advantages  over  equations  (33)  through 
(35).  The  system  can  support  steady  solutions;  this  formal  conclu¬ 
sion  is  consistent  with  Miller's  void  observations,  which  showed  a 
steady  void  distortion  in  the  nutating  system.  This  system  can  be 
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easily  linearized  and  the  conditions  for  linearization  are  obvious 
(see  Section  5)#  Equation  (39)  displays  an  important  forcing  term, 
2rflr  =  -2tor  cos  $ ,  which  cannot  be  incorporated  into  the  reduced 
pressure.  In  dimensional  form,  this  term  is  proportional  to 
uifi  sin  e.  By  comparison  with  equation  (23),  it  becomes  clear  that 
this  term  produces  the  deviation.  In  fact,  the  system  equations 
(38)  through  (41)  have  a  trivial  solution  if,  and  only  if, 

Re  •  flr  =  0.  In  the  following  sections,  we  continue  the  analysis 
based  on  the  equations  and  notation  introduced  for  the  nutating 
system . 


4.2.2  The  Reduced  Pressure  p*. 

The  relation  between  the  reduced  pressure  pn  and  the 
pressure  p  is  given  by 


Pn  =  P  -  2  Un  x  I 


1  2 


P  ‘  2  1 


(42) 


where  r,  is  the  distance  from  the  nutation  axis 
d 


From  equation  (37)  we  obtain 


P'  *  F  *  2 


(  1  +  )  r2  +  X2 

z  d 


(43) 


Using  equation  (36),  we  find  for  point  r,  <p,  z 


P  =  P  “  2 


1  +  T  (  1  -  0 


1/212 


J  r^+(Torsin4>)^ 


+  (T0z)^  +  2  T  ^  cr  (  l  -  a  ^  rz  cos  <£ 


( 


(44) 


It  is  straightforward  to  show  that  p’  is  identical  with  the 
reduced  pressure  p^  in  the  body-fixed  system  for  the  same  point. 
The  difference  p  -  p'  has  to  be  taken  into  account  in  calculating 
forces  normal  to  the  cylinder  walls  and  for  the  moments 
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L4  ’ 


1 

■ 


U 


* 

► 

w 


perpendicular  to  the  z-axis.  No  contribution  is  made  to  the  despin 
moment  about  the  z-axis,  that  arises  from  only  tangential  shear 
stresses  on  the  walls. 

5.  LINEARIZED  EQUATIONS 

The  momentum  equation  [see  equation  (39)]  is  nonlinear  due 

to  the  convective  terms  contained  in  D' .  There  is  little  hope  of 

finding  solutions  of  the  full  system  equations  [see  equations  (38) 

through  (41)]  except  approximations  from  computational  or 

perturbation  methods.  Perturbation  methods  often  enhance  insight 

into  the  structure  of  the  problem  and,  therefore,  we  consider  in 

this  section  the  lowest-order  approximation  for  small  deviations 

v  ,  v.  ,  v  ,  p'  from  rigid-body  motion,  small  enough  for  neglecting 
r  o  z 

the  quadratic  nonlinear  terms. 

There  are  two  circumstances  that  lead  to  small  deviations. 
The  first  case  of  weak  forcing  through  ftr  was  mentioned  above.  The 
second  case  occurs  at  a  small  Reynolds  number,  that  is,  large  vis¬ 
cous  damping.  The  motion  is  then  essentially  governed  by  the  terms 
multiplied  by  1/Re  in  equation  (39).  We  will  consider  this  case  in 
Section  5.3 

5 . 1  The  Case  of  Weak  Forcing . 


for 


According  to  equations  (39)  and  (36),  weak  forcing  occurs 


to  <<  1  ( 45 ) 

This  condition  can  be  satisfied  by  either  small  t  (nutation  rate 

small  in  comparison  with  the  spin  rate)  or  small  a  (small  nutation 

angles),  or  a  mixture  of  both.  It  is  not  necessary  to  restrict 

2 1 

and  o  separately  as  was  done  in  some  previous  studies. 
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With  equation  (45)  satisfied  and  velocity  components  of 
order  0 ( T <f> )  #  equations  (39)  can  be  reduced  to 


3  v  3  v 

inr  +  2(1  +  ^lz)  v<<> 


,  ,  r  v  „  3v. i 

=  _  i£l  +  1_  D»  v _ £  _  2 _ £.] 

3r  Re  [D  r  r2  r2  34>  | 


(46a) 


ii. . 

3t  34> 


2d  +  n  )  vr 


=  -  -  +  ^-  D" ' 

r  3$  Re 


V  3v  T 

_  Jl  +  2_ 

2  2  3<|> 

r  r  J 


(46b) 


3v  3v 

_z  .  _ z 

"St  34> 


-  2rnr  -  js  D-vz 


(460 


Continuity  equation  (38)  and  boundary  conditions  equation  (41) 
remain  unchanged.  This  system  of  equations  is  linear,  but  of  high 
order  due  to  the  presence  of  viscous  terms.  It  is  the  key  to 
understanding  the  moments  exerted  at  the  cylinder  for  arbitrary 
Reynolds  number  and  weak  forcing.  The  steady  solutions  at  small 
values  of  Re  are  of  special  interest.  A  first  attempt  to  find 
these  solutions  in  closed  form  led  to  yet  unresolved  problems  with 
the  boundary  conditions  for  a  finite-length  cylinder  We  return  to 
this  point  in  Section  5.3.  Only  for  small  viscosity,  that  is,  for 
the  limit  Re  »  ,  the  system  equations  (38  ),  (  46),  and  (41)  allow 
for  relatively  simple  solutions. 


5. 2  The  Inviscid  Limit. 


For  Re  00  and  equation  (45),  continuity  equation  and 
momentum  equations  take  the  following  form: 
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(47a) 


I  !_ 

r  3r 


(rvr>  +  r  1? 


9z 


3<J> 


2(1  +  “z1  v4>  =  "  if” 

2<*+«z> 


(47b) 


(47c) 


3v 


z 


at 


+ 


-W- 


-  2rfl. 


( 47d ) 


This  system  is  similar  to  that  considered  by  Stewartson  and 
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Roberts,  who  neglected  flz.  Using  an  ansatz  of  the  form 


vr  =  fr(r'z>  sin  (  4>  +  kt),  V(f  »  f^  (r,z)  cos  (  $  +  kt) 

v  =  f  (r,z)  sin  (  $  +  kt),  p*  =  g(r,z)  cob  (  <t>  +  kt)  (48) 

z  z 

one  obtains 


(1  +  k)  fr  -  2(1  +  Slz)  f$  «  -  If 

-d  +  k)  f$  +  2d  +  nz)  fr  =  7  g 


(49a) 

(49b) 


(1  ♦  k)  fz 


2r 


to  cos  0 _ 

cos  TT+XtT 


The  first  two  equations  provide 

_  J-u  -  2U  +  y  a] 

r  ~ (1  +k)2-4(l+n)^ 

z 


(  49c) 


(50a) 
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(50b) 


-2(1  +a2)  %  -  a  +  n  f 

(i  +  k>2  -  4(1  +  n,)2 


Continuity  equation  (47a)  requires 
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* ,  j 


+  I  la 

3r2  r  9r 


a_  a  +  *>2  -  id  ♦  nz)2  329 
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(1  +  k) 


3z‘ 


which  can  be  written  as  a  Laplace  equation  for  g,  by  suitably 
transforming  z.  A  particular  steady  solution  of  equation  (51)  for 
k  *  0  can  be  found  in  the  form 


90  ■  Arz  +  Br,  p'Q  =  gQ  cos  $ 
Consequently 


(52) 


v 


r 


0 


( Az  +  B)  sin  $ 


_  ( Az  B )  COS  <P 

- I~r^r 


v  *  (2xa-  A)  r  sin  i 
z0 


(53) 


This  solution  accounts  for  the  steady  forcing  term  in  equation 

(49c).  Additional  unsteady  solutions  with  k  /  0  can  be  found  from 

equations  (49),  (50),  and  (51)  with  the  forcing  term  in  equation 

(49c)  dropped.  These  unsteady  solutions  represent  the  modes  of 

inertial  oscillations.  Under  more  restrictive  conditions,  these 

modes  were  studied  by  Stewartson  for  a  cylindrical  container,  and 

11  21 

by  Stewartson  and  Roberts  for  a  spheroidal  container.  '  The 
spheroidal  container  is  an  interesting  special  case  where  only  one 
mode  of  oscillation  of  the  liquid  is  induced  by  nutation. 


... 


M 
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For  the  spheroidal  container  with  aspect  ratio  X 


r2  +  -^  =  1  (54) 

.  X2  ■  • 

considered  by  Stewartson  and  Roberts,  the  boundary  conditions  can 
be  written  in  the  form 


rv 


r 


0 


(55) 


It  follows  that  B  *  0  in  equation  (53),  but 

A  =  2to(l  +  2nz)/(l  +  -  X2)  (56) 

The  unsteady  solutions  are  free  oscillations  of  arbitrary  amplitude 
C,  which  is  fixed  by  initial  conditions  and  frequency: 


k 
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1  +  2flz  -  X 


1  +  X 


(57) 


If  container  shape  X  and  nutation  rate  and  angle  are  such  that 
k  y  0,  the  solution  is 


2%z 

v  =  - s  +  Cz  sin  (  $  +  kt) 

r  1  +  2HZ  -  X 


-2nrz 


1  +  2Q  - 

Li 


— t  +  Cz  cos  (  <f>  +  kt ) 
X^ 


v 

z 


-2  X2fi*r  , 

- -  -  CX^r  sin  (  4>  +  kt) 

1  +  2  n  -  x 
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it 
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» 


a. 


p' 


-2(1  +  2 flz)  nrrz 

1  ♦  2°*  -  *2 


\2c  +  n  ) 

+  2C  - 5 —  rz  cos  ( $>  +  kt) 

1  +  r 


(58) 


For  a  sphere,  X  »  1,  the  frequency  from  equation  (57)  is  k  =  ft_. 
This  result  is  different  from  Stewartso  i  and  Roberts,  who  found 
resonance,  k  =  0  for  X  =  1.  This  result  is  reproduced  by  letting 
ft_  0  in  equations  (57)  and  (58). 

Here,  resonance  occurs  in  a  slightly  nonspherical 
2 

container  with  X  =  1  +  ft  .  For  k  -  0,  we  obtain  a  solution  in  the 

Z 

form 

Az 

vr  ■  r-is;  sin  ♦  -  rm;  “ 

%  -  rna;  c°?  ♦  *  2t 

(1  +  2ft  )  ft 

v2  =  -Ar  sin  +  +  - m -  rt 

z 

p’  =»  Azr  cos  $  +  —  "-r ft-  j^ftr  -  (1  +  2ftg)  ft^tj  (59) 

where  A  is  an  arbitrary  constant  governed  by  the  initial  condi¬ 
tions.  The  physical  interpretation  of  the  solutions  (see  equations 
(58)  and  (59)]  is  very  similar  to  that  given  by  Stewartson  and 
Roberts . 


Disregarding  the  constants  in  these  solutions,  which  are 
as  yet  undetermined,  we  find  that  velocity  components  and  the 
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pressure  p'  are  proportional  to  to,  or  of  order  O(tct).  We  also 
note  the  Simple  periodicity  of  the  solution  in  the  azimuthal 
coordinate  <J>. 


We  have  not  yet  found  analytical  solutions  of  the  viscous 
equations  (46)  subject  to  equations  (38)  and  (41).  With  the 
solutions  of  the  previous  section  and  those  of  Stewartson  for  the 
inviscid  case,  however,  we  can  obtain  some  qualitative  information 
on  viscous  solutions.  ^  The  inertial  oscillations  will  suffer 
increased  damping  as  the  Reynolds  number  decreases.  The  resonant 
peaks  in  the  side  moments  versus  T  will  broaden  and  finally  perish. 
At  a  sufficiently  low  Reynolds  number,  the  equations  support  a 
steady  solution.  The  dependence  of  this  solution  on  <p  will  still 
be  simple.  However,  the  simultaneous  appearance  of  velocity  com¬ 
ponents  and  their  derivatives  with  respect  to  <p  in  equation  (46), 

2 

for  example,  3vr/3$  and  -vr/(Re  •  r  )  in  equation  (46a)  indicates  a 
phase  shift  of  the  velocity  field  with  respect  to  the  forcing. 

This  shift  increases  as  Re  decreases. 


Exploiting  the  simple  periodicity  in  <J>,  the  linearized 
equations  can  be  reduced  to  partial  differential  equations  for 
functions  of  r  and  z.  A  separable  solution  can  be  found  only  for 
an  infinite  cylinder,  X  ->■  <*>.  We  have  not  pursued  this  solution  for 
two  reasons.  First,  it  is  questionable  whether  this  solution  is 
relevant  to  the  problem  with  X  =  0(1).  Second,  the  solution  of  the 
resulting  ordinary  differential  equation  plus  boundary  conditions 
requires  a  major  computational  effort.  For  finite  aspect  ratios, 
the  main  difficulty  for  the  analytical  work  is  to  satisfy  the 
boundary  conditions  at  side  walls  and  end  walls  simultaneously. 
Similar  problems  were  found  by  Gerber,  Sedney,  and  Bartos.^  They 

employed  the  boundary-layer  approximation  in  order  to  satisfy  the 
no-slip  condition  at  the  end  walls.  Therefore,  the  range  of 
applicability  of  their  method  is  comparable  to  Murphy's  method, 
which  uses  t lie  boundary-layer  approximation  at  end  walls  and  side 


walls.  ^  These  mechods  are  not  applicable  at  low  Reynolds  numbers 
where  viscous  boundary-layers  thicken  and  finally  merge, 
overwhelming  the  inviscid  core.  For  this  viscous  range,  further 
analytical  work  needs  to  be  carried  out.  Alternatively,  the 
partial  differential  equations  for  the  functions  of  r  and  z  can  be 
solved  by  computational  methods.  The  problem  of  solving  the 
linearized  equations  appears  as  a  basically  two-dimensional  special 
case  of  the  full  problem  based  on  the  nonlinear  equations  in 
Section  4.2 

5.4  Scaling  Aspects. 

From  equations  (38),  (41),  and  (46),  we  have  made  various 
attempts  to  obtain  modified  systems  of  equations  that  reveal  the 
functional  dependence  on  combinations  of  the  dimensionless 
parameters.  We  have  rot  yet  been  successful  in  identifying  single 
(or  paired)  parameter  combinations  that  govern  the  solution.  The 
dependence  on  X  concealed  in  the  boundary  conditions  can  be 
introduced  into  the  differential  equations  by  stretching  the 
z-axis,  z  =  Xz .  In  order  to  keep  the  continuity  equation  free  from 
parameters,  the  axial  velocity  can  be  replaced  by  v  *  Xv  .  The 
parameters  t  and  o  introduced  by  the  forcing  term  in  equation  (46c) 
can  be  incorporated  in  the  velocity  scale.  The  form  of  X^  in 
equations  (53)  or  (58)  suggests  rescaling  the  velocity  components 
and  pressure  p '  by  a  factor Xto,  which  would  replace  the  reference 
velocity  ua  by  ftc  sin  0.  We  note,  however,  that  no  change  of  the 
Reynolds  number  would  occur  in  the  linear  equations  for  the 
deviation  from  rigid-body  rotation.  Additional  information  on 
proper  scaling  is  expected  from  pursuing  analytical  solutions  of 
these  equations. 

5 . 5  Estimates  on  the  Velocity . 

In  order  to  obtain  some  guidance  on  the  validity  of  the 
linearized  equations,  we  consider  Miller's  experimental  data,  in 
particular  his  figure  13.^  With  0  =  20°  we  obtain  a  =  0.342,  while 


T  varies  between  0.125  and  0.25.  Therefore,  0.043  <  tc  <  0.086, 
which  can  be  considered  small.  It  is  questionable,  however, 
whether  the  resulting  deviations  from  rigid-body  motion  are  small 
enough  for  linearization.  With  X  =  4.29,  the  inviscid  forced 
solution  scales  with  Xto  £  0.367,  which  may  be  too  large  to  permit 
use  of  equation  (46).  on  the  other  hand,  viscous  damping  will 
reduce  the  deviation. 

The  measured  despin  moments  allow  a  rough  estimate  of  the 
velocity  gradients  at  the  cylinder  walls.  These  moments  originate 
from  the  shear  stresses  at  the  side  wall,  and  t2^,  at  the  end 

walls.  We  assume  that  tZ(^  and  xr^  are  of  the  same  order  and  can  be 
replaced  Dy  an  average  value  x.  The  despin  moment  can  then  be 
written  as 


The  term  1/3  represents  the  contribution  from  the  end  walls  and  can 
be  dropped  in  comparison  with  X  =  4.29,  so  that  equation  (60) 
contains  only  the  average  of  x^  for  r  =  a.  This  stress  is  defined 
as 


r<J> 


■  vp 


r=a 


"  3v  3v,  v,‘ 

I 

r  3<p  3r  r 


r=a 


With  vr  =  v<j,  *  0  at  the  side  wall,  we  obtain 


(61) 


T  = 

X 

=  vp 

■SvT 

r<p 

t  =a 

3r 

r  =  a 


(62) 


Exploiting  equation  (60),  the  average  gradient  of  at  the  side 
wall  can  be  expressed  as 
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v/a  = 
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(63) 
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With  M  =  0.7  N  •  m  for  v  =  1  m  /sec,  and  p  =  1400  kg/m  ,  the  aver¬ 
age  velocity  gradient  is  v/a  =  0.04/sec,  providing  an  average 
velocity  of  the  order  v  =  0.002  m/sec.  The  dimensionless  value 
with  respect  to  the  angular  velocity  of  the  rigid-body  rotation 
( ui  =  3000  rpm)  is  v/(wa)  =  0.00013,  very  small  indeed.  For 
v  =  0.01  m  /sec,  the  values  are  v/a  =  8.5/sec,  v  =  0.5  m/sec, 
v/(wa)  *  0.027,  still  sufficiently  small  to  allow  using  the  linear 
equations.  Therefore,  the  linear  equations  seem  to  cover  the  full 
range  where  large  despin  moments  were  observed.  For  the  smaller 
viscosity  v  =  10  m  /sec,  the  velocity  gradient  assumes  a  large 
value,  v/a  =  254/sec.  Formally,  we  obtain  v  =  15.3  m/sec.  This 
number  may  be  misleading,  however,  since  the  growth  of  the  velocity 
occurs  only  across  the  thickness  o  <  a  of  the  viscous  boundary 
layer.  The  estimates  from  equation  (63)  provide  some  guidance  on 
the  magnitude  of  the  velocity  only  in  the  fully  viscous  regime,  say 
for  Reynolds  numbers  Re  <  )000. 


6.  REMARKS  ON  MORE  GENERAL  CASES 

From  the  foregoing  discussion,  it  appears  that  the  viscous 
fluid  motion  in  the  spinning  and  nutating  cylinder  can  be  studied 
on  the  basis  of  the  linearized  equations  of  Section  5.  However, 
there  are  a  number  of  disconnected  topics  that  require  using  other 
sets  of  equations.  Concerning  the  gyroscope  experiments  and  the 
still  unexplained  changes  in  the  yaw  angle  growth  rates,  the 
solutions  of  the  inviscid  equations  should  be  reconsidered  without 
resorting  to  weak  forcing  (equation  (45)].  Beyond  the  inertial 
modes  studied  by  Stewartson,  which  have  period  2 tt  in  <J> ,  other 
classes  of  modes  exist  with  period  2iT/n,  n  >  1,  which  may  be 
excited  through  the  generation  of  harmonics  at  larger  yaw  angles.^ 
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Other  cases  of  interest  are  the  viscous  fluid  motions  at 
smaller  spin  rates,  where  equation  (45)  may  not  be  satisfied. 
Whenever  equation  (45)  is  invalid,  the  products  of  flr  and  with 
velocity  components  in  equation  (39)  must  be  retained.  Even  at  a 
small  Reynolds  number,  these  products  generate  harmonics  in  <J>,  and 
the  flow  field  assumes  the  form  of  a  Fourier  series  in  $ . 

Only  the  linearized  equations  of  Section  5  bear  the  promise 
of  analytical  or,  at  least,  semianalytical  soxutions.  Strong 
forcing  or  large  deviations  from  rigid-body  motion  can  only  be 
treated  by  numerical  methods.  There  is  as  yet,  little  experience 
in  computing  internal  flows  in  simply  connected  rotating 
containers.  The  three  problem  areas  that  need  to  be  overcome  are 

•  No-slip  condition  at  the  end  walls 

•  Three-dimensionality  of  the  flow  field 

•  Failure  of  most  methods  at  the  axis  r=0. 

As  mentioned  above,  use  of  the  linearized  equations  reduces 
the  problem  to  two  dimensions,  r  and  z.  We  have  studied  the 
applicability  of  three  classes  ot  methods — finite  difference, 
polynomials,  and  finite  element. 

Finite-difference  methods  seem  to  be  least  promising  due  to 
the  difficulty  of  achieving  sufficient  accuracy  in  presence  of  the 
axis  at  r  =  0,  The  simple  geometry  of  the  cylinder  favors  using 
spectral  methods,  with  a  Fourier  expansion  in  4>  and  a  Chebyshev 
expansion  in  z.  Use  of  Chebyshev  polynomials  in  r  was  considered, 
but  abandoned  for  slow  convergence  of  the  series.  Legendre 
polynomials  or  Jacobi  polynomials  in  r  appear  as  the  better  choice, 
but  additional  numerical  studies  are  necessary  for  evaluating  the 
convergence  properties.  Whereas  polynomials  are  usually  chosen  for 
computational  ease,  expansions  in  Bessel  functions  may  work  best 
for  the  present  problem.  It  requires,  however,  major  efforts  to 


develop  this  yet-unexplored  branch  of  spectral  expansions.  The 
third  method  that  seems  most  straightforward  to  apply  to  the 
present  problem  is  a  finite-element  method.  No  difficulties  are 
encountered  on  the  axis.  Even  for  the  three-dimensional  problem,  a 
relatively  small  number  of  elements  may  provide  sufficient  accuracy 
of  the  moments.  The  accurate  calculation  of  the  flow  field  with  a 
large  number  of  elements  is  a  straightforward  extension  at 
increased  computational  expense. 
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